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The pseudogap phase of the underdoped cuprates is the host to a variety of novel electronic 
phenomenon. An example is the dc-resistivity which shows metallic behaviour in the ab-plane, 
while the c-axis response is insulating. We apply a model, originally formulated by Yang, Rice, 
and Zhang, to study the resistivity in the pseudogap phase. This model is able to reproduce the 
qualitative features of the resistivity, including the systematic deviations from linear behaviour for 
the in-plane conductivity, and the insulating behaviour along the c-axis. We compare this to the 
predictions of the arc model and find similar qualitative behaviour. We find that the most important 
element in understanding the resistivity is the reconstruction of the Fermi surface, which puts strong 
restrictions on the number of quasiparticles allowed to participate in dc-transport. 



The psuedogap is the precursor phase from which 
high-T c superconductivity condenses in the underdoped 
cuprates. Understanding the pseudogap phase is believed 
to be central for deducing the mechanism that underpins 
high-T c superconductivity To study the low lying 
excitations in this phase, it is essential to understand 
the character of its Fermi surface (FS). Photoemission 
(ARPES) in the psuedogap phase observes segments of 
FS along the Brillouin Zone diagonals. These segments 
were originally interpreted as disconnected segments of 
FS and became known as Fermi arcs Q. Recently, 
they have been resolved as pockets with small spectral 
weight along the aniferromagnetic Brillouin zone (AFBZ) 
boundary [3J. Evidence for pockets is compounded by 
the observation of quantum oscillations, which requires a 
closed FS 0] ■ The nature of the FS reconstruction in the 
pseudogap phase has been the subject of much debate 
and has lead to a number of theoretical proposals. 

One class of theories of the pseudogap phase involve 
competing order parameters whose fluctuations serve to 
reconstruct the FS into the observed pockets. These 
models involve broken translational symmetry, such as 
charge or spin density waves. Other theories of the pseu- 
dogap begin from the idea of preformed Cooper pairs 
that lack off-diagonal long range order Indeed, su- 
perconducting fluctuations have been seen to exist above 
T c in the form a disordered vortex liquid [6J. However, 
these fluctuations fail to persist to the temperature scale 
associated with pseudogap behaviour Q. An alterna- 
tive approach is that the pseudogap arises naturally as 
one dopes a Mott insulator, which then leads to both 
insulating and strongly correlated electronic behaviour 
Q. This was the approach that Yang, Rice, and Zhang 
(YRZ) adopted when they put forward their model of the 
underdoped cuprates in 2006 

The YRZ model consists of an ansatz for the coherent 
piece of the single particle Greens function in the pseu- 
dogap phase. It is based on results for Anderson's res- 
onating valence bond (RVB) spin liquid 13, 
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FIG. 1. (Color online) In plane (solid blue, in units e~ 2 d) and 
c-axis (dashed orange, in units e~ 2 t^ 2 d~ 1 ) resistivities as a 
function of temperature in the YRZ model for dopings x = 
0.05, 0.10, 0.15 and 0.2. The solid turquoise lines are fits to 
the linear T region. The c-axis resistivity shows insulating like 
behaviour, and the in-plane resistivity shows a systematic de- 
viation from linearity induced by the pseudogap energy scale. 
The insets show the Fermi surfaces (shown for one quarter of 
the BZ, AFBZ boundary indicated by the dashed line) given 
by the YRZ Greens function. 



description, as the Mott insulating state is approached a 
gap opens on the AFBZ boundary. This gap is an RVB 
spin gap and appears as an energy scale separate from 
superconductivity. 

Since its debut, the YRZ model has been shown to cap- 
ture many properties of the cuprates that were considered 
anomalous. In particular, it has been shown to give good 
agreement with Raman Spectra 12, "HI. ARPES (3, 14|, 
specific heat 1151 . penetration depth [16j . and tunnelling 
spectroscopy [17] . Perhaps most remarkably is that the 
YRZ model captures all of this behaviour with doping as 
the only free parameter. 

A related model, the arc model, has also been shown 
to capture the qualitatively correct behaviour in the un- 
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derdoped cuprates (see for example Refs [13|, In 
this model, a gap is placed on the antinodal portion of the 
metallic fermi surface, leaving Fermi arc segments. The 
length of the remaining FS serves the same role doping 
does in the YRZ model. Thus, it seems that the qual- 
itative behaviour of many of these experiments can be 
captured by any model that gives the disappearance of 
FS along the antinodal direction. It is crucial to under- 
stand in what circumstances the YRZ and arc models 
provide qualitatively different pictures of the electronic 
response in the pseudogap phase. 

In this letter, we compute the in-plane (ab), and out- 
of-plane (c-axis) resistivity using both the YRZ model 
and the arc model for the pseudogap phase of the under- 
doped cuprates. This calculation is of an essentially dif- 
ferent character than the previous studies, as it requires 
us to add ingredients not present in the initial formula- 
tion of YRZ. Specifically, the resistivity relies on the scat- 
tering rate, and the extension to the c-axis needs some 
assumption to be made about the interlayer tunnelling 
matrix element. After adding these ingredients we show 
that the YRZ and arc models agree qualitatively with 
dc-transport experiments. Although these two models 
have different underlying physics, they give qualitatively 
similar predictions for the resistivity in the underdoped 
cuprates. 

YRZ Model — The YRZ ansatz for the coherent part 
of the electron Greens function for a doping x is 
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where gt = 2a;/ (1 + x) is a Gutzwiller rcnormaliza- 
tion factor, £(k) = £°(k) — At'(x) cos(k x ) cos(fc y ) — 
2t" '(x)(cos(2k x ) + cos(2k y )) — [i p is the tight binding dis- 
persion on a square lattice out to third nearest neigh- 
bours, and £°(k) = —2t(x)(cos(k x ) + cos(k y )) is the first 
nearest neighbour term. /i p is a chemical potential deter- 
mined from the Luttinger sum rule. The RVB gap, Apg, 
plays the role of the pseudogap and has d-wave symme- 
try, i.e., A PG = A (l - x/x c )(cos{k x ) - cos(k y )). Wc 
use values for all parameters in the YRZ model as they 
appear in Rcf. [9[. We work in units where H = fee = 1> 
and all energies are measured in units of to (see Ref. [§]). 

Equation (fT]) would be the Greens function for a su- 
perconductor if £°(k) = £(k). This difference causes the 
gap to open on the surface defined by £°(k), which in 
our case is the AFBZ. This gap naturally reconstructs 
the large FS into pockets and can give rise to interesting 
electronic response. Given this Greens function we can 
compute the conductivity, and hence, the resistivity. 

The dc-conductivity is given by the Kubo formula 
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uated at the one-loop level 

IIy(w) = e 2 T K(k)G(k, v n ) Vj (k)G(k,u + v n 



(3) 



The current-current correlation function, 1X^ (0;), is eval- 



where vi is the z th component of the velocity and the 
v n are Matsubara frequencies. For the ab-conductivity 
we use v x — d£,(k)/dk x . For the c-axis conductivity, we 
replace v 2 with i^(k)d 2 , where d is the interlayer dis- 
tance and t±(k) = t± (cos(k x ) — cos(k y )) 2 is the inter- 
layer tunnelling matrix element. This matrix element 
was first suggested by Anderson It was used in the 
nodal liquid model of the cuprates, which reproduced 
the insulating-like behaviour associated with the pseudo- 
gap [1 91 ] - The final component that we need to compute 
the conductivity is the scattering rate, T, which broad- 
ens the spectral densities. To make contact with experi- 
ments we follow the suggestion of reference [2(| and take 
a temperature dependent scattering rate. For our calcu- 
lations we include a small residual scattering and take 
T = 0.01 + 2tt\T, with A = 0.3. These choices leave us 
with doping as our only tuneable parameter. 

Figure [1] shows the in-plane and out-of-plane resistivity 
as a function of temperature for a few values of doping. 
The qualitative agreement with experiments is good for 
both the in-plane and out-of-plane results (see references 
2oL 2l|). At high temperatures the in-plane resistivity 
is linear. As the pseudogap opens (moving down in dop- 
ing), there is a progressively larger deviation from the 
linear behaviour. This deviation from linearity is a loss 
of metallically associated with the shrinking of the FS. 
The c-axis resistivity shows strong insulating behaviour 
at low values of doping, and becomes increasingly metal- 
lic as optimal doping is approached, just as in the exper- 
iments. 

To understand these results, we found it instructive to 
introduce a joint density of states (JDOS). The JDOS can 
be thought of as a reference function, and its computation 
is the same as Eq. ([3]), but with the velocity operators 
omitted. This function is shown for the optimally doped 
and underdoped cases in figure [5] for the c-axis. The 
only difference in the JDOS between the ab- and c-axis 
conductivity is the factor (cos(fc K ) — cos(fc ?/ )) 2 from the 
tunnelling matrix element. This factor suppresses the 
spectral weight near zero frequency. 

At optimal doping the JDOS contains only a single 
peak, and the dependence on the resistivity is controlled 
solely by thermal factors. At low temperatures, when 
the thermal factors are more sharply peaked, we get a 
higher conductivity and hence a lower resistivity. As we 
increase the temperature the thermal factors broaden, we 
get a smaller conductivity (larger resistivity) . The under- 
doped case is more interesting. In the underdoped case 
the opening of the pseudogap causes the peak to split 
into two pieces, separated by an energy on the order of 
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FIG. 2. (Color online) Joint density of states involved in the 
conductivity calculation for x = 0.05 (highly underdoped) 
and x = 0.20 (optimal doping) at low and high temperatures. 
The red dashed curves are the thermal factors for the low 
temperature case and peak height ~ 25. The blue dashed 
curves are the the thermal factors for the high T case. As the 
JDOS becomes more flat (psedogap fingerprint blurred out) 
we return to a more metallic-like behaviour. 



the gap scale. The thermal factors behave as before, but 
there is no spectral weight remaining in the low T case. 
This gapping of the region near uj = naturally gives 
insulating behaviour. As the temperature is increased 
the JDOS broadens due to increased scattering, and our 
resistivity falls back down accordingly. Within the YRZ 
model the resistivity is determined by this interplay be- 
tween thermal factors and the structure of the density of 
states coming from the pseudogap. 

Arc Model — We now turn our attention to the arc 
model. In the arc model, we start with a large FS as in 
the optimally doped case. By hand we place a gap, A, 
on part of the FS starting from the antinodal direction. 
The amount gapped out is paramaterized by an angle, 
8 C , which plays the role of doping. There are then two 
contributions to the conductivity: a 'free electron' part 
(fare) f rom the remaining FS, and an 'interband' part 
(<ta) from the gapped region of the FS. We then com- 
pute the conductivity for a scattering rate T for both the 
in-plane, and c-axis. We obtain (units as in Figure 3) 
(afc) and = 
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FIG. 3. (Color online) The in-plane resistivity (units of 
SmA > ^ * s ^ ne density of states at the Fermi energy, and Vf 

is the Fermi velocity.) and c-axis resistivity (units of e 2^ t i ) 
in the arc model for dimensionless pairing strengths n — 0, 
0.25, 0.5, 0.75, and 1. 
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where A = A/%/ A 2 + T 2 , i^xIA 2 ) is the incomplete ellip- 
tic integral of the first kind, and E(x\X 2 ) is the incom- 
plete elliptic integral of the second kind. 

Unlike in YRZ, in the arc model the arc length is in- 
dependent of the magnitude of the pseduogap. To make 
contact with the experiments and YRZ we chose to trade 
8 C for a new dimensionless variable k — A/Ao, where Ao 
is the magnitude of the pesudogap at doping. In prin- 
ciple the function 9 c (k) is a strictly increasing function of 
k which could be obtained by fitting the arc model to ex- 
perimental data, if one desired. To highlight the essential 
features, and for simplicity, we chose 8 C — kit/ 4. The re- 
sulting curves are shown in figure [3l Since the scattering 
rate T oc T, the x-axis may be read as temperature. 

Although the precise details of the curves differ from 
both the YRZ model and the experiments, the broad 
features are still captured by this simple model. The in- 
plane resistance is linear at high T, and shows systemati- 
cally larger deviations as the paring strength is increased. 
The case k = 1 is distinct from the other cases, and 
corresponds to the universal limit first discussed by Lee 
[22| . The c-axis resistivity shows progressively insulat- 
ing behaviour as we proceed further into the underdoped 
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FIG. 4. (Color online) The c-axis resistivity with the weight 
given by the remaining FS as a tuneable parameter, a. Specif- 
ically, a enters as acr^ c + . This figure shows how the FS 
contributes to the c-axis conductivity. The FS still gives an 
important contribution, even though the tunnelling matrix 
element to gaps out the nodal excitations. 



regime. This is due to the loss of FS along the antinodal 
direction. The c-axis tunnelling matrix element is then 
responsible for ridding of the remaining (nodal) states, 
which leads to an insulating like behaviour. 

The results of the arc model can be understood simply 
as a two oscillator model: one Drude at w = 0, and a 
displaced oscillator at u = A. The dc-conductivity is the 
value from these two contributions at uj = 0. As the scat- 
tering rate is increased these two contributions broaden. 
At first, the broadening causes a decrease in the conduc- 
tivity as the Drude peak comes down. This continues 
until the displaced peak at the gap energy leaks into the 
lj = region, leading to a gain in conductivity. The con- 
ductivity then tends to decrease as the peaks broaden 
further. This explains the weak maxima seen in both the 
ab- and c-axis resistivity (Figure [3]). To show the two 
contributions explicitly we weighted the FS contribution 
to the conductivity by a parameter < a < 1 (Figure 
[4]) . As we weaken the Drude component of the conductiv- 
ity, the sample becomes increasingly more resistive. This 
depletion of the remaining quasiparticles is precisely the 
same physics behind the resistive behaviour coming from 
the c-axis tunnelling matrix element. 

In conclusion, we have found that the YRZ model of 
the pseudogap state with a linear in temperature scatter- 
ing rate naturally reproduces the resistivity seen in this 
phase. A simpler model, the arc model, was also shown 
to possess qualitatively correct features of the resistivity. 
The physics in the two models which gives rise to the re- 
sistivity are very different. In the YRZ model, it was the 
interplay between thermal factors and a density of states 



which was responsible for the behaviour of the resistiv- 
ity. By contrast, the density of states in the arc model 
is constant and temperature only entered through the 
scattering rate. However, owing to the two distinct con- 
tributions to the conductivity, the arc model still showed 
qualitatively correct features. 

Two elements in the YRZ calculation are vital for the 
agreement with experimental resistivity. First is the re- 
construction of FS as a function of doping. This re- 
construction reduces the number of quasiparticles avail- 
able for dc-transport as the Mott insulating phase is ap- 
proached. The second element is the tunnelling matrix 
element. This matrix element effectively removes the re- 
maining nodal quasiparticles and leaves an insulating re- 
sponse. The fact that FS reconstruction is so important 
allows us to understand why the arc model can do well, 
despite lacking all of the microscopic elements of the YRZ 
model. Indeed, one can think of the arc model as a phe- 
nomenological version of YRZ, where one is willing to 
ignore the mechanisms behind the reconstruction of the 
FS. Nevertheless, the arc model can still be used to help 
deduce the underlying physics, and provides simple ana- 
lytical results which are of great use. 
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